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the relationship of red and photographic 

INFRARED SPECTRAL DATA TO GRAIN YIELD 
VARIATION WITHIN A WINTER WHEAT FIELD 


ABSTRACT 

Two-band hand-held radiometer data from a winter wheat field, collected on 21 dates during 
the spring growing season, were correlated with within field final grain yield. Significant linear rela- 
tionships were found between various combinations of the red and photographic infrared radiance 
data collected and the grain yield. The spectral data explained -64% of the within field grain yield 
variation. This variation in grain yield could not be explained using meteorological data as these were 
similar for all areas of the wheat field. Most importantly, data collected early in the spring were 
highly correlated with grain yield; a five-week time window existed from stem elongation through 
an theses in which the spectral data were most highly correlated with grain yield; and manifestations 
of wheat canopy water stress were readily apparent in the spectral data. 
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THE RELATIONSHIP OF RED AND PHOTOGRAPHIC 
INFRARED SPECTRAL DATA TO GRAIN YIELD 
VARIATION WITHIN A WINTER WHEAT FIELD 

REPRODUCIBILITY OF THE 

BACKGROUND ORIGINAL PAGE IS POOR 

Several approaches for remote sensing of winter wheat yield have been proposed, The Large 
Area Crop Inventory Experiment (LACIE) involving NASA, NOAA, and the USDA, was perhaps 
the most visible and the best known of various winter wheat yield prediction approaches. Simula- 
tion or regression models were used to predict wheat, yield based upon climatic conditions (MacDonald 
and Hall, 1977). 

The EarthSat Corp, (1976) lias proposed a spring wheat yield prediction method. This approach 
differs from the LACIE approach in that NOAA meteorological satellite data were used to estimate 
grid cell precipitation which in turn drove the yield determination. Critical variables for this approach 
were the calculation of soil moisture and grid cell weather from the meteorological satellite data and 
the determination of crop phenology ( LACIE, 1978). 

Idso et al. (1977a,b) have proposed a method of winter wheat yield prediction using the stress 
degree day concept, That is, the final yield of a crop was hypothesized to be linearly related to the 
accumulated stress degree days over some critical period. This technique was developed in Phoenix, 
Arizona under irrigated conditions and is currently being evaluated in dryland win ter wheat growing 
areas. 

Hammond (1975) and Moram and Williams (1975) have discussed monitoring wheat production 
with satellite remotely sensed data. Harlan and Liu (1975) and Colwell et al. (1977) have evaluated 
the use of Landsat data for inferring direct winter wheat grain yield predictions. Spectral data were 
used to estimate yield related variables such as stand density or leaf area index, Because grain yield 
is usually highly correlated to stand density or leaf area index, a direct estimate of final grain yield 
was usually possible. Colwell et al. (1977) concluded that Landsat (i.e,, spectral) data explained a 


■considerable amount of yield variation which was not explained by meteorological data. In addition, 
Landsat derived yield predictions were as highly correlated with individual field yields as were esti- 
mates made using traditional sampling techniques, even, in some cases, if the Landsat data were 
collected several weeks before the field samples. 

Heilman et at. (1977) reported a technique for estimating winter wheat grain yields using Land- 
sat derived estimates of leaf area index coupled with an evapotranspiration model, More recen tly, 
Wiegand et al. (1979) discussed leaf area index estimates for winter wheat made from Landsat data 
and the implications for evapotranspiration and crop simulation modeling using these data, 

Accurate inputs for crop canopy leaf area index are a necessity for successful evapotranspiration 
and crop simulation modeling, To achieve accurate leaf area index estimates from Landsat data, 
“ground-truth” sampling must occur which adequately samples the variability of the actual field leaf 
area index for enough Landsat pixels ( — 0.45 ha) to satisfy basic sampling theory requireme nts. This 
has often proven difficult to achieve. Heilman et al. (1977), Kanemasu et al, (1977), and Wiegand 
et al. (1979) all used three ground samples of 91 cm by 2 row widths (~40 cm) to establish each field’s 
(>40 ha) leaf area index. Tims ~1 m 2 was used to determine the leaf area index of a large field. This 
could be a significant source of error for the evapotranspiration and crop yield modeling approaches 
(Wiegand et al., .1979). “Ground-truth” sampling must adequately account for the in situ population 
variability by inclusion of enough samples to allow for a reasonable estimate of the actual field leaf 
area index. 

DESCRIPTION OF RESEARCH UNDERTAKEN 

As part of an ongoing research program into future satellite sensor development including selec- 
tion and evaluation of spectral bands, radiometric resolution, frequency of coverage, orbit selection, 
and other considerations for vegeta tional applications, wo have been collecting hand-held radiometer 
data from a variety of agricultural crops in a NASA/GSFC-USDA/BARC cooperative research program. 
Previous research on alfalfa, corn, and soybeans had demonstrated that red and photographic infrared 
two-band radiometer data were highly correlated with various properties of plant canopies (Tucker 
et al, 1979a,b,c). Therefore, we decided to evaluate the applicability of these data for yield predictions 
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on wheat. A companion paper reports on the use of these same data for monitoring total dry mat- 
ter accumulation (Tucker et al,, 1979d). 

Previous work with the green leaf or photosynthetically active biomass has suggested that this 
physiological entity integrated the various biotic and abiotic effects present in the plant canopies 
(Tucker et ah, 1 973; Colwell et ah, ) 977; Tucker ct ah, 1 979a, b). Conditions which adversely af- 


fected plant growth and development resulted in a reduction in the photosynthetically active biomass. 
Because the photosynthetically active biomass or green leaf area is one of the basic system variables 
in primary production, monitoring this system variable throughout the growing season should enable 
inferences to be made regarding total dry matter accumulation and grain yield, Tins lias been pro- 
posed as the Leaf Area Duration (LAD) concept (Colwell et ah, 1 977 ; Richardson et ah, 1 979), 

However, the LAD concept needs to be modified in that the photosynthetically active biomass 
or leaf area is actually the interaction between the green LAI and the chlorophyll concentration. 
Expressed in other words, the photosynthetically active biomass can be defined as the interaction 
between inter- and intra-leaf scattering and chlorophyll absorption which occurs predominately in the 
green leaves of the plant canopy in question. 

Red and photographic infrared spectral data have been demonstrated by many workers to be 
highly correlated with the photosynthetically active biomass of several cover types (reviewed in 
Tucker, 1979), The red spectral data are highly correlated witli the in vivo chlorophyll concentra- 
tion, whereas the photographic infrared data are highly correlated with LAI. Thus, various linear 
combinations of these two adjacent spectral regions are highly related to the photosynthetically active 


biomass (Tucker, 1979). 
EXPERIMENTAL PROCEDURE 
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Our experiment was conducted in a 1 .2-ha soft red winter wheat ( Triticum aestivum L.) field 
at the Beltsville Agricultural Research Center, Beltsville, Maryland. The field was plowed, disked, 
and planted with the cultivar ‘Arthur* on October 6, 1977 at a rate of 107.6 kg/ha. A conventional 
grain driil with 1 7.8-cm row spacing was used for seeding. Before seeding the field was limed on the 
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basis of soil test recommendations and fertilizer was applied at a rate of 33.3 kg N, 53.8 kg l>, and 
53.8 kg K/lia. The following spring (early March 1978) the crop was topdressed with 20.4 kg N/lia. 

Twenty 2- x 3-m plots were selected during the winter dormant period In the wheat field. Sub- 
sequently, for each plot, four pairs of red (0.65*0.70 pm) and photographic infrared (0,775-0.825 pm) 


spectral radiance measurements were obtained using a hand-held digital radiometer (Pearson et ah, 
1976) on 21 dates between March 21, 1978 (Julian date 80) and June 23. 1978 (Julian date 174). 


The intervals between dates ranged from I to 9 days; however, the average interval was 4.7 days and 
the median interval was a tie between four and five days (Table I), 

The red and photographic infrared spectral radiance data were used to form the ir/red ratio and 
the normalized difference (ND) of Rouse et al. (1973) and Dee ring et ah (1975) where: 

ND = (IR - RED)/(IR + RUD) (I) 

The four pairs of the spectral measurements per plot were averaged to account for the spatial 
variability present in each plot. All spectral data were collected plus or minus 90 minutes of local 
solar noon, measured normal to the ground surface at a height of ~1 m above ihe plant canopy under 
sunny skies (Table 1 ), 


Throughout the growing season average plant height, percentage cover estimates, and phenol og- 
ical development notes were recorded for the field area (Table 2). The crop reached harvest maturity 
in Late June 1978. On June 28, 1978 (Julian date 1 79) a 0.9- x 3, 0-m swath was cut with a small 
sickle bar mower from the center of each plot and the grain thrashed with a small plot grain thrasher. 
The grain yields were oven dried at 60°C for 72 hours and were subsequently adjusted to 14% mois- 
ture and expressed in g/m 2 . 

A regression approach was taken in which the average red radiance, ir radiances, ir/red ratio, 
and ND were correlated With the grain yield for each of the 21 data collection dates. In addition, 
the ir/red ratio and the ND as a function of Julian date and growing degree days were integrated for 
four different time intervals and correlated with grain yield. 
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Table 2 

Agronomic Data Pertaining to Average Plant Heights, Estimated Percentage 
Canopy Cover, and Crop Growth Stages at 1 1 Selected Dates for the 
20 Winter Wheat Plots from 1978 


Calendar 

Date 

r 

Julian 

Date 

Plant 

Height 

(cm) 

Percentage 

Cover 



, 

i- 


04/24/78 

114 

35.0 

54 

05/01/78 

121 

45.2 

56 

05/11/78 

131 

70.8 

66 

: 05/19/78 

139 

90.8 

64 

05/25/78 

145 : 

112.0 

68 

06/01/78 

152 : 

112.5 

61 

06/06/78 

157 ; 

114,8 

63 

■■ 06/14/78 

i 

165 

115.5 

64 

. 06/20/78 

17! ! 

111. 8 

51 

06/23/78 

174 j 

108.5 

51 

; 06/27/78 

178 | 

104.0 j 

51 


Growth Stages 


Numerical 
(after Xndoks 
etfll., 1974) 

Descriptive 

34 

stem elongation. 4th 
node detectable 

35 

stem elongation, 5th 
node detectable 

44 

booting, boots just visible 

58 

inflorescence emerges 

64 

an thesis, half-way 

73 

early milk 

85 

soft dough 

85 

soft dough 

87 

hard dough 

89 

hard dough 

92 I 

ready for harvest 


RESULTS AND DISCUSSION 

Examples of the radiance data were plotted against Julian date (Figure 1), The red radiance 
vs, Julian date curve showed the influence of increasing chlorophyll absorption to ~ Julian date 


139 followed by a gradual decrease in chlorophyll absorption as a result of the onset of senescence 
(Figure la). The photographic infrared radiance vs. Julian date curve showed a gradual increase 
with time to a maximum at '-Julian date 139 followed by a decrease caused by senescence. This 
resulted from the direct relationship of the photographic infrared radiance to the foliar density of 


the plant canopy (Figure lb). The ir/ red radiance ratio and the normalized difference both exhibited 
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40 

JULIAN DATE 


RIvMvOUUCIHIUTy OF TIIE (B) 
original PAGE IS poor 



120 140 

JULIAN DATE 



120 140 

JULIAN DATE 


Hgure 1, Rod radiance, photographic infrared radiance, ir/red radiance ratio, and the 
normalized difloicncc plotted against Julian date for one ol 20 wheat plots sampled. 
(A) Red radiance, (B) Photographic infrared radiance, (C) Ir/rod radiance ratio, and 
(D) Normalized dillerencc (ND). Note how the ir/red radiance ratio and the normal- 
ized difference effectively compensate for the variability in the radiance data. 






1 similar trends with respect to Ju'. i date (I gures Ic and Id, respectively). Both spectral variables 

showed the same five component trends with respect to Julian date They both increased as the 
*pring portion of the growing season began until -Julian date 102 where their rate of increase de- 
clined: front Julian date - 102 to - 1 1 2 both showed rapid rates of increase; from Julian date —112 
to - 1 23 there w as a decrease for both spectral variables; from Julian date - I 23 to - 130 there was 
another increase in both spectral variables; and from Julian date - 139-175 there was the senescence 
caused decrease in both the it red radiance ratio and the normalu.J difference. 

I flicsv same trends were observed for the majority of our experimental plots and were an « sample 
ot a manifestation ot transient wneat canopy water stress (I mure 2). The ir/rvd radiance ratio and the 
normalized difference both exhibited a decrease in their rate ot increase by Julian date - 102 This 
decrease was ended by a period of rainfall on Julian date 102 which, coupled with warmer tempera- 
tures. resulted in a rapid rate of increase in the ir/ted radiance ratio and normalized difference until 
Julian date - 1 1 2 where a period of mild drought stress again occurred, reaching maximum drought 
stress severity on - Julian date 1 23 The later episode of drought stress w as ended by a period of week- 
long rains which prevented any field data collection lor the interval of Julian dates 125-130 The 
plant canopy subsequently recovered from the previous water stress condition and reached maximum 
values for the ir red radiance ratio and the normalized difference on Julian date 139 which corre- 
sponded to lull spike emergence i Table 2 1 From t his date on, the progression of wheat canopy senes- 
ccncc resulted in .» reduction ot the ir red radiance ratio and the normalized difference (Figure 2). The 
spectral manifestation of drought stress was consistent with previous research results (Thompson and 
Wchmancn. 1979 Tucker, et al 1979c). Figure 2 also demonstrated the dynamic nature of plant cano- 
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pics.tnd how rapidly they respond to and recover from water stress conditions. We interpreted the mech- 
anism for this spectral manifestation of plant canopy water stress to be largely due to a reduction in the 
leaf chlorophyll density . Iliis followed from Figure la where there were decreases in red absorption 
(i.e.. inen .ises in radiance) lor Julian dates !()2 and 1 23 which corresponded to the episodes of water 
stress. I lie photographic infrared radiance, by contrast, did not show an\ effect which was attributed 


* - ~- 
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Figure 2. The 
2- x 3-m plots ] 
rainfall. Note I 


JH/RED RADIANCE RATIO 



to the episodes of water stress (Figure 1 b), This implied that the LAI remained fairly constant while 
the leaf chlorophyll density was temporarily reduced, either through photooxidation, enzymatic, or 
some other reducing mechanism (Tucker et ah, 1973; Tucker et al., 1 979c), 

The next phase of the analysis regressed each of the 4 spectral variables against the final grain 
yield for each of the 21 data collection dates. Our principal attention was focussed on the ir/red 
ratio and the normalized difference as these linear combinations adjust for irradiational variability 
(Figure 3). We observed lower correlations for the first four dates, higher correlations for dates 5 
and 6, and a marked decrease in correlation with grain yield for date 7 which was attributed to very 
windy conditions (Tables 1 and 3), Sampling dates 8, 9, and 10 were all highly correlated to grain 
yield. There was a slight decrease for date 1 1, followed by high correlations for sampling date 13 
after which the correlations decreased as senescence progressed (Figure 4). 

We observed a 40-day period (from Julian date 112-1 52) during which our spectral data were 
more significantly correlated to the final grain yield than earlier or later (Figure 4; Table 3). How- 
ever, the regression relationships were not constant during this period, suggesting that a regional 
application of these relationships was doubtful (Table 4), An alternative to tills was attempted in 
the form of integrating the spectral data in terms of Julian date or accumulated temperature units 
(growing degree days or equivalent). This was a remote sensing application of the leaf area duration 
concept. 

The ir/red radiance ratio and tile normalized difference were integrated for four periods during 
the growing season: Julian dates 80-112; 1 12-152;152-174; and 80-174. We observed that the plateau 
portion (Julian dates 1 12-1 52) of the growth curves, corresponding to the maximum green leaf bio- 
mass present, was the most highly correlated with final grain yield (Figure 5). The normalized differ- 
ence spectral data from this period explained 66% of the within field variability in final grain yield while 
the normalized difference spectral data from the entire growing season (i.e., Julian dates 80-174) ex- 
plained 64%of the same grain yield variability (Figures 5 & 6). Spectral data alone thus explained ap- 
proximately two-thirds of the variability in within field grain yield. This could not be explained by mete- 
orological models as tire micro-climatic effects were largely identical for our 1.2-ha field. The earlier 
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Correlation Coefficients for the Four Radiance Variables and Final Grain Yield 
for Each of the 21 Days where Spectral Data were Collected from 20 2 x 3 Plots 
in 1978, Refer also to Figures 3 and 4. 

Sampling 

Sequence 

Julian 

Date 

Ret! 

Radiance 

Ir 

Radiance 

Ir/ Rod 
Radiance 
Ratio 

Normalized 

Difference 

1 

80 

0.14 

0.51* 

0.27 

0.27 

2 

89 

-0,33 

0.48* 

0.41 

0.43 

3 

92 

-0.22 

0.42 

0,31 

0,30 

4 

95 

-0.30 

0,43 

0.40 

0.42 

5 

97 

-0.7 1** 

0.78** 

0.62** 

0.75** 

6 

102 

-0.7 1*"' 

0,80** 

0.61** 

0.75** 

7 

104 

-0.51* 

0.59** 

0.55* 

0,57** 

8 

112 

-0,75** 

0,74** 

0.69** 

0.78** 

9 

1 18 

-0.82** 

0.66** 

0.73** 

0.82** 

10 

121 

-0.82** 

0.54* 

0.76** 

0.82** 

1 1 

123 

-0.78** 

0.64** 

0,68** 

0.79** 

12 

13) 

-0,76 

0.78** 

0,75** 

0.78** 

13 

139 

-0,78** 

0,83** 

CO 

oo 

o 

0.81** 

14 

146 

-0,73** 

0.76** 

0.77** 

0.76** 

15 

152 

-0.78 

0.69** 

0.78** 

0,78** 

16 

157 

-0.73** 

0.61** 

0.71** 

0.73** 

17 

161 

-0.55** 

0.68** 

0,66** 

0.66** 

18 

165 

-0,51* 

0.49* 

0.56** 

0.55** 

19 

166 

-0.57** 

0.43 

0.61** 

0.58** 

20 

170 

0.19 

0,51* 

0.32 

0.30 

21 

174 

0.46* 

-0.14 

-0.48* 

-0.42 


Indicates significance at the 0.05 level of probability 
Indicates significance at the 0.01 level of probability 
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Tabic 4 

Regression Derived I stimales for d () and p’| lor Julian Dales I I 2-1*2 
Using the Ir Red Radiance Ratio and the Normalized DifTerenee to Predict 
the I'inal Grain Yield. The I quation used was of the I orm Grain Yield - 
do 4 h> | (Spectral Variable* 


Sampling 

Julian 

It Red Radiance i 

latio 

Nut mali/ed Difference 
i 1 

Period 

Date 

00 




0. 


8 

112 



0.48 


365 0 

0,61 

9 

IIH 

230 0 


0.53 


350 3 


10 

121 

224.4 


0.57 


327.0 

0.68 

II 

123 



0.46 


321 1 


i: 

131 

207.7 




413 


13 

139 

173.1 


0.69 


511 


14 

146 

170.8 

21.5 

0.5‘> 


475. ~ 


15 

152 

148.4 

29.0 

O.o 1 


534 A) 

0.60 



1 
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(B) 

Figure 5 The r 2 values resulting from the (A) ir/red radiance ratio ~ Julian 
date and (B) normalized difference - Julian date integrations for four differ- 
ent time periods: Julian dates SO-1 ] 2; 1 1 2-1 52: 1 52-1 74; and 80-1 74 f Re- 
fer also to Figure 6). ' 
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integration period of Julian date 80-1 12 only explained 0,48# of the variability in final grain yield 
Figure 5). 

Any successful application of these relationships for regional remote sensing prediction of 
winter wheat grain yield must be able to account for within region differences in winter wheat crop 
phenology. Horn Table 4 we showed that the regression derived equation coefficients varied for 
the 40-day period (Julian dates 1 1 2-152) where highly significant correlations between the spectral 
data and final gram yield were reported (Tables 3 and 4). Thus, slight crop phenology differences 
within a given region would preclude the effective use of one m two satellite images obtained during 
this 5- to 6-week per iod. hither accurate crop calendar information must be available or the leaf 
area duration technique must bg employed. 

Crop calendar information is frequently lacking in extra-country situations. A possible remotely 
sensed .substitute for this could be the combination of spectral data and accumulated temperature 
units (figure 7). rim iollows from the basic relationship of temperature to biological activity. 

When combined with red and photographic infrared spectral data, representing tire plant canopy 
vigor or potential lor growth, this may have between region and/or between year comparison utility. 



li EPPOE UCIBI I TTY OF THE 
ORIGINAL PAGE IS POOR 


1920. 

GROWING DEGREE DAYS 


3200. 


Figure 7, 
and 


A possible spectral crop calendar using tile normalized difference 
accumulated temperature units (i.e., growing degree days). 
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In the event that a combination of spectral anti accumulated temperature unit data would be 
unsatisfactory, frequent satellite coverage would be required to apply the leaf area duration concept 
for growing season integration ol spectral data, As shown in Figure 2, frequent satellite coverage is 
needed to record the rapid onset of plant canopy stress conditions. It is easy to visualize the situa- 
tion in Figure 2 if there were only 4, 5, or 6 data points instead of the 21 we collected. Obviously, 
entire stress periods might not he detected if only a few observations wore available. For yield con- 
siderations, the occurrence oi the stress condition is of paramount importance, Ifa period of water 
sti css pectus during heading or during the grain tilling period, the reduction of the grain yield is much 
greater than if this same stress condition occurs at some other time. 

Spectral data are capable of providing large-area information highly related to crop condition 
or vigor, Whether these data are used in conjunction with other remotely sensed information, by 
simulation modelers, or by themselves, it is apparent that remotely sensed red and photographic 
infrared spectral data can supply valuable crop condition information. 
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CONCLUSIONS 

' hrCqUently COllcc!ed spcctral lllUa wcrc sh0W1 ' to be highly related to crop condition. Two peri- 
ods ol water stress were readily apparent, 

2. A 40-day period extol from ft. beginning of stem elongation through anthe* during which ft, 
speemd data were highly correlated with final grain yield and explained at least 60* of the variability. 

3. The integrated spectral data explained 64* of the variability In within field final grain yield. This 

could not be explained using meteorological models as the mlcro-elhnatc conditions were very similar 
lor our experimental wheat field, 

4. I-arge-seale extra-count, y application for our findings would require „ spectral crop calendar. We 
propose a possible spectral-acciinil, luted temperature unit crop calendar, If this is not feasible, wc 
suggest frequent satellite observations and application of the leaf area duration concept. 

5. The combination of spectral data and other environmental or agronomic dad, will in, prove the 
predictive Utility of yield forecasting. While spcctral data are useful by themselves, they are more 
useful in combination with other types of data. 
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